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Abstract 

We show that rotondike excitations are thermally induced in a two-dimensional dilute Bose gas 
as a consequence of the strong phase fluctuations in two dimensions. At low momentum, the 
roton-like excitations lead for small enough temperatures to an anomalous phonon spectrum with 
a temperature dependent exponent reminiscent of the Kosterlitz-Thouless transition. Despite the 
anomalous form of the energy spectrum, it is shown that the corresponding effective theory of 
vortices describes the usual Kosterlitz-Thouless transition. The possible existence of an anomalous 
normal state in a small temperature interval is also discussed. 
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I. INTRODUCTION 



Superfluidity in interacting Bose systems has been a fascinating research topic for several 
decades. After the achievement of Bose-Einstein condensation in dilute atomic gases,— many 
of the remarkable properties of superfluid 4 He were also observed in these weakly interacting 
systems, such as vortices and Josephson oscillations. 2 In order to understand their superfluid 
properties, it is necessary to clarify the role of the elementary excitations in determining 
these properties. According to a celebrated criterion due to Landau, a Bose gas without 
interactions cannot be a superfluid, since its spectrum e p = p 2 /2m makes it indifferent 
to Galilei transformations. A superfluid, on the other hand, is resistant to a slowdown of 
the molecules due to the loss of Galilei invariance. A Bose system with the relativistic 
looking excitation spectrum E p = cp has this property, thus being a superfluid by Landau's 
argument. Thus, phonon-like excitation spectra are an essential part of a superfluid. 

In superfluid helium, which is a strongly interacting Bose system, the interatomic po- 
tential has an attractive short-range part over a distance do, the average interatomic dis- 
tance. This sets the scale of a further important set of excitations. Scattering experi- 
ments of neutrons show that around a momentum p pa 1/ao the spectrum behaves like 
E p = A + (p — po) 2 /2m*. These elementary excitations are called rotons. 

Feynmani was the first to recognize the importance of the rotons for a superfluid. In 
modern language, Feynman's theory describes rotons as a result of large quantum phase 
fluctuations at low temperatures. These produce small vortex loops of size ~ clq. At higher 
temperatures where thermal fluctuations take over, the vortex loops combine to larger ther- 
mally excited vortex loops which profit from the high configurational entropy of line-like 
excitations. At the critical A-point, these loops become infinitely long and destroy the order 
of the superfluid.^ 

In two dimensions, phase fluctuations are so strong that they destroy the long-range order 
at any temperature.^ As noted first by Kosterlitz and Thouless (KT)£ there still exists a 
phase transition driven by phase fluctuations. At low temperature, a film of superfluid 4 He 
contains vortex- antivortex pairs bound by Coulomb attraction, whose unbinding causes the 
KT transition. One of the most important predictions of this theory is the universal jump 
to zero of the superfluid stiffness at the critical temperature.— This destroys the excitations 
of energy E = cp, and thus the superfluidity. 
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The increased relevance of phase fluctuations in two dimensions suggests that Feynman's 
rotons are more abundant than in three dimensions. We want to argue that this is indeed 
the case: in spite of the weak interactions, a dilute Bose gas possesses roton-like excitations 
which are the precursors of the high-temperature vortices that lead to a KT phase transition. 
This is quite remarkable since, contrary to superfluid helium, the bare interaction does not 
contribute with an attractive part to fix the size of po- 

The plan of the paper is as follows. In Sect. II we briefly review the known results of 
the t-matrix formalism in d dimensions for dilute Bose gases. Sect. Ill contains most of 
the main results of the paper. There we apply the so called dielectric formalism^ to the 
two-dimensional dilute Bose gas. We will essentially work out a RPA approximation in two 
dimensions. Since the dielectric formalism is well known in three-dimensional applications, 
we will not discuss its derivation here, referring the reader instead to the literature. The 
focus will be in the application of the method to two dimensions. From our calculations a 
new, temperature dependent, excitation spectrum emerges, namely, 



We will show that the above spectrum allows for thermally excited roton-like excitations and 
that its low momentum behavior exhibits an anomalous power behavior with a temperature 
dependent exponent t]o{T). Remarkably, in the two-dimensional dilute Bose gas at finite 
temperature, the RPA analysis will lead to closed form analytic results including roton-like 
excitations. In Sect. IV we discuss the effect of phase fluctuations in the system. Our 
analysis in this Section will allow us to derive the approximate critical temperature of the 
system as r)o(T c ) = 1/4. The actual critical temperature is determined in the usual way 
following the Kosterlitz-Thouless vortex unbinding mechanism; 6 It will also be shown that 
a crossover temperature T* exists, above which our anomalous spectrum becomes unstable. 
Sect. V concludes the paper. The whole discussion is supplemented with two Appendices 
containing calculational details. 

II. THE EFFECTIVE INTERACTION AT ZERO TEMPERATURE 

In general the effective interaction for any dimension d < 4 depends on a momentum 
scale p associated to the energy of the scattered particles. At zero temperature the effective 




(1) 
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interaction in d dimensions reads 

ATi d l 2 a d - 2 /m 

9 ^ ~ 2 2 - d T(l - d/2){pa) d ~ 2 + T(d/2 - 1) ' ^' 

where m is the reduced mass of the scattered particles and a is the s-wave scattering length. 
The above result is essentially exact at zero temperature and corresponds to a geometric sum 
of ladder diagrams which gives the only nonvanishing contribution to the vertex function.^ 
From Eq. (j2J) we derive the exact renormalization group (RG) /5-function for the dimen- 
sionless coupling g(p) = mp d ~ 2 g(p) 

d n-d/2) „ 2 



9 ~"~ 2 d + 1 7T d / 2 



(3) 



For 2 < d < 4 the only fixed point is g* = 0. On the other hand, for d < 2 a nontrivial fixed 
point located at 

2 d+1 TT d / 2 

9 * ~ ~dT{-d/2) [ ) 

exists. Thus, d = 2 is the upper critical dimension for the dilute Bose gas. Writing d = 2 — e 
and expanding Eq. (j3J) for small e, we obtain 

g* ps 2ne. (5) 

In the dimension interval 2 < d < 4 we can easily take the limit pa — > in Eq. (J2J to obtain 

4:TT d / 2 a d - 2 

90 ~ r(d/2-l)m" (6) 
For d — 3 the above equation reproduces the familiar formula go = 4na/m. For d — 2, 

however, we cannot set pa = 0. In two dimensions Eq. (0) becomes 

/-x 2n/m 
92DiP) = ln2- 7 -ln(pa)' (?) 
where 7 is the Euler constant. Eq. (JJ|) is in agreement with previous work . n i 12 

The first theories for the two-dimensional dilute Bose gas were developed by Popov 13 
and Schick.— More recently, an improved version of Popov's theory was presented by Stoof 
and collaborators . 15 ! 16 In Popov's approach,— the bare interaction go<5( x ) i s replaced by an 
effective interaction g5(x) determined by a t-matrix, leading in d dimensions with d € [2, 4) 
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Arr d / 2 a d - 2 /m 

9d ~ 2 2 - d T(l - d/2){na d ) d - 2 + T{d/2 - 1) ' (8) 



where n is the density and a the s-wave scattering length. The above interaction corresponds 
precisely to the one given in Eq. (J2J), where we have set pa = na d . For d-»2we obtain the 
two-dimensional coupling constant of the Popov-Schick theor y 13 i 14 

,. 2n/m 
^"i im ^ = ~ ln(eW/2) ' (9) 

where 7 is the Euler-Mascheroni constant. The logarithm in the denominator implies that 

the effective repulsion decreases only very slowly with decreasing density.^ Fisher and 

Hohenbergii have shown within the Popov-Schick theory that the dilute limit na d <C 1 

of the d > 2 theory must be replaced for d = 2 by lnln(l/na 2 ) ^> 1. 



III. THE DIELECTRIC FORMALISM IN d = 2 

The t-matrix result (j2J) incorporates, via a Lippmann-Schwinger integral equation, the 
sum of all ladder diagrams. In this Section we take into account the sum of bubble diagrams 
of the plasmon type,— which are nonvanishing at finite temperature. This corresponds to 
the random phase approximation (RPA), which sums geometrically the diagrams shown in 
Fig. 1. This approximation has often been applied in d = 3 dimensions.^ In RPA the vertex 
function containing the effects of both resummations is given explicitly in d dimensions by 

r(iw lP ) = - 9 4~. — -, (io) 

1 - gdRfa, p) 

where the polarization bubble is given by 
~ 1 \ f d d q 

U(iuj, p) = ~22 J J^y G °( iu} + iuJn > p + q) G o(^., q), (11) 

with G (iuj,p) = — e p ). The chemical potential is canceled by the Hartree 

contribution.-- 7 The denominator of (jllj) determines the regular part e r (iu),p) of the di- 
electric function e(iu),p) renormalizing the density correlation function, i.e., 

^^-ik^ (12) 

where 

e(iu, p) = 1 - g d [Tl (iuj, p) + U(iu, p)], (13) 

with 

n ° (!w ' p)= («^r (14) 
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FIG. 1: Feynman diagram representation of the vertex function Eq. (|1U|), The dashed line rep- 
resents the bare interaction. The vertex function is obtained as a geometric series of polarization 
bubbles. 

contributing to the singular part of the dielectric function.— The regular contribution, being 
given by a particle-hole diagram, vanishes at zero temperature. We will treat U(iui, p) in 
the so called classical approximation, where the Bose-Einstein distribution function tib{x) = 
\j[eP x — 1) is replaced by ub{x) ~ T/x. Thus, the classical approximation is valid in the 
limit p\x <C 1, where At = (2mT)~ 1 ^ 2 is the thermal wavelength. Therefore, we obtain in 
the limit p\x <C 1 the result 

U(iu, p) = -A d m 2 Tp d - 4 [ e -( d " 2 )/ 2 (1 - tuj/e^ + (1 + iuj/e^] , (15) 

with A* = 2 2 ~ d n- d / 2 T(d/2-l)T(3-d). A derivation of Eq. (HSJ) using Feynman parameters 
is given in the Appendix. For d = 3, this reduces to the well-known result^ 

d=3 _.Tm 2 , fiu + e p 



U(zu;,p) a ^-z— -In . (16) 

2-np \iuj — e p J 

In d — 2 dimensions Eq. (fTB^) has a pole associated with a logarithmic short distance 
divergence. We will remove this divergence using the s-wave scattering length as short- 
distance cutoff. This leads to 

~ r sd=2 2mT e p \n(pa) 

U(iaj,p) = —— -. 17 

7r {luoy - e 2 

The spectrum of elementary excitations is obtained from the poles of Xppi^, p) . 9 i 19 i 20 At zero 
temperature, only IT (iu;,p) contributes, and we recover the usual Bogoliubov spectrum. 




FIG. 2: The RPA corrected excitation spectrum of Popov-Schick theory exhibiting a roton-like 
minimum. The parameters in a system of units such that ks = fr = 1 are m = 0.5 cm -2 sec, 
n = 0.01 cm -2 , T = 1.2 sec -1 , and a = 2 cm. 

The excitation spectrum is obtained from the pole of the density correlation function, 
which corresponds to the vanishing of the dielectric constant (fH?j) . In this way we obtain 
the excitation spectrum announced in the Introduction: 



E p = \\ei + 2g 2 ne p 



Tm 

1 In (pa) 

7m 



(18) 



The above energy spectrum possesses a thermally induced roton-like minimum. The excita- 
tion spectrum is shown in Fig. 2 for a suitable set of parameters. 

The approximate position of the roton minimum can be determined in the following way 
using the Landau criterion for superfluidity. According to Landau criterion, the critical 
velocity above which excitations appear in the fluid is given by the minimal value of the 
ratio E p /p. This minimum value corresponds to the point p = po for which 

dEp _ Ep^ 
dp p 

On the other hand, the energy spectrum (fTHj) satisfies the equation: 
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A -^ = T + wX E ^^r)- (20) 



The point po for which Eqs. and (|2~U|) coincide is determined by the vanishing of the 
second term on the right-hand side of Eq. (J2(Jj) . i.e., 



Po(T) = J- = mJ^. (21) 

Qt \ TV 



The above value of p corresponds approximately to the position of the roton-like minimum 
of the excitation spectrum (|I8jl. The T-dependent length scale ax nearly replaces the a of 
bulk 4 He. 

The consistency of the calculation leading to the spectrum (J18)) can be checked by com- 
puting the spectrum also from the pole of the anomalous propagator. It is well known 
that the pole of the propagator should be the same as the one from the density correlation 
function,«^*22, although most approximations fail to fulfill this requirement. After analytically 
continuing to real frequencies, we obtain that the pole of the propagator is given by the 
solution of the equatior 8,!) 



UJ 



2 2 



e 2 p -2ne p T(uj,p) = 0. (22) 



The above equation is a generalization of the Bogoliubov result for the excitation spectrum. 
Indeed, it corresponds to an improvement of the Bogoliubov result where the coupling con- 
stant #2 is replaced by the vertex function (fTUj) . 8 Solution of Eq. (j22j) with the polarization 
bubble (JT7J) gives precisely the energy spectrum (JTSJ). 

The excitation spectrum pSjl has another interesting property in the low-momentum 
regime, which for mg^j^ < 1 is defined by 

p> « p2( T ) = ^]f£9l < 1 = 2mT 4?m . (23) 
71 A T 

In the above low-momentum regime the logarithm gives rise to an anomalous power behavior 



' 92 n a -Tm/2Trn^l-Tm/2irn 

m 

The exponent of p in <^^ can be rewritten as 



x m 



" (T)5l -S = 1 "» ,n (25) 

where p s (0) = mn is the superfluid mass density at zero temperature, and r] (T) = 
Tm? /27rp s (0). Interestingly, precisely the same exponent arises in the spin wave treatment 
of classical phase fluctuations in the study of the KT transition. We will revisit this analysis 
in the next Section, where it will also be shown that despite the anomalous scaling of the 
spectrum, the KT transition occurs as usual at higher temperatures, in agreement with the 
discussion of Ref . Q. This is to be expected, since the KT transition actually occurs through 
a vortex-antivortex unbinding mechanism. 
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The elementary excitations described by Eq. (}2~4")) are stable only for rj (T) < 1, i.e., 
for T < 2iip s (0)/m 2 = T*. Note that the low- momentum inequality (J2~3|) already requires 
T < T». But this condition can be softened to just T < in which case Eq. (J24*j) is still 
approximately valid up to logarithmic corrections. The role of the temperature T* will be 
discussed later. 

It is instructive to calculate the low-momentum contribution to the superfluid density 
following from (|24jl . Since vortices are not included in the above calculation, the result will 
show only the anomalous phonon contribution to the so called background superfluid density, 
p s o- The calculation is based on the Landau prescription,— according to which the superfluid 
density is p^{T) = p — pJ^(T), where p£o(T) is the anomalous phonon contribution to the 
normal fluid density. From the Landau formula the normal background fluid density is given 
by 

(3 f d d p p 2 e^ 
Pn0 ~ dj (27r) rf (e^-l) 2 ' 1 ] 

Now we set d = 2 and insert the anomalous phonon spectrum (}2"4"|) in Eq. (|26|) to obtain 



n ph (T) w - 1 ' V 



4 ~ g(T) 
a(T) 



C 



4 ~ or(T) 
a(T) 



(27) 



where c = (g2n/m) 1 ^ 2 a~ rio( - T \ By expanding this in powers of T, we obtain to leading order 
the usual phonon contribution, which in two dimensions is proportional to T 3 . In deriving 
Eq. (|27|) we have assumed the usual hydrodynamic limit where the upper cutoff — here 
Po(T) — is taken to be infinity. 



IV. PHASE FLUCTUATIONS 



By integrating out small density fluctuations in the hydrodynamic limit, we obtain the 
following effective action for the phase fluctuations: 

(28) 

where the effective Hamiltonian contains a local and a nonlocal interaction between the 
superfluid velocities v s (t, r) = W(r, r)/m: 

tj rjlocal i rj-non— local /on\ 

H cS = H c{i + H cS , (29) 



s, 



off 



dr 



/ d 2 r{d T 9) 2 + H cS 

2 92 J 
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where 



H%* = ™ I dVv s 2 (r,r), (30) 

^non-local = 1 ^ ^ /" _ ^ ^ p ) . ^ ^ (3^ 

being a bilocal mass density, which is obtained from the Fourier transform of Ai(p) = 
nm(pa)~ 2ri0 ^ T \ In the folowing discussion we will show that at high temperatures the non- 
local part contributes only in a small temperature interval above T c . Below T c the usual 
effective Hamiltonian for the phase fluctuations given by Eq. (jHO)) dominates the critical 
behavior and the KT transition obtains. The arguments to be described below consider the 
scaling behavior of the spin wave theory for /j°° n ^ local anc l the field theory for the vortices, 
which consists of a generalized sine-Gordon theory. 

Let us consider first the case without the RPA correction, i.e., in the absence of the non- 
local effective Hamiltonian. This just corresponds to the usual phonon spectrum. In such a 
situation the effective Hamiltonian is given simply by Eq. The spin wave analysis is 

in this case well known £2i However, it is useful to review it here in order to compare with 
the non-local spin wave regime. 

At higher temperatures we can neglect the higher Matsubara modes so that d T 9 = and 
the effective action becomes 

n i „ — >-)2 



and 



with 



ScS = 2riTJ d r[w(r)] ' (33) 

and the problem is essentially a classical one. Let us recall the computation of the correlation 
function (ip(r)i(j*(r')) in this regime and in the absence of vortices (spin- wave theory) In 
such a case we just have to compute the correlation between the phases: 

( e i[0(r)-0(r')]} = _>_ y D^-J^'fe[V9(r")] 2 +»J(r'')0(r'')} ) (34) 

where J(r") = 5 2 (r" — r) — 5 2 (r" — r'). The Gaussian integral is straightforward and yields 

( e Wr)-«(rO]) = exp |^[G(r - r') - G(0)]} , (35) 



where 



J2 tpr 

CM = I V- <36) 
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In order to evaluate the above Green function we introduce the regularized Green function 

f d 2 p e ip r 

GM2{r) = ] J^f^TW^ (37) 

where M is a regularizing mass to be sent to zero at the end of the calculations. Evaluating 
the integral explicitly, we obtain 

G M <r) = ^K (Mr), (38) 
where K (y) is the modified Bessel function of the second kind. On the other hand, we have 

G M2 (0) = ^lnfAy (39) 



2vr \M, 

where A = n/a is the ultraviolet cutoff. Now we can safely take the limit M — > 0: 

G(r) - G(0) = lim [Gm 2 (r) - G A/2 (0)] = In (-) + const. (40) 



Therefore, 



<e*>-<«> ~ r _ 1 (41) 



r — r 



Note that precisely the exponent t]q(T) that we have defined in Eq. (|23j) arises in the 
above equation. However, the mechanism that generates the anomalous behavior in the 
above classical spin-wave theory is completely different from the quantum case discussed 
in the previous Section. Indeed, there the anomalous scaling of the spectrum arises due to 
interaction effects, while in the above calculation it follows from the analytic properties of 
the Green function of a Gaussian classical theory in two dimensions. 

Let us now study the classical problem associated to the effective Hamiltonian at large 
distances. In this case (pHj) dominates, since the corresponding power of p is smaller than 
two. Once more we neglect the higher Matsubara modes to obtain the effective action for 
the classical problem as 

S c s = ^ [ d 2 r [ d\'M(r - r') v s (r) • v s (r'), (42) 



2T „ 

where v s (r) = v s (0, r). The correlation function between the phases has again the form 
(J35J) . except that the Green function is now given by 

f d 2 r> p ipr 

G « - 1 WfTMfy (43) 
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We have 



g2r)o(T) r2ir poo ^ipr cos 



G(r) = Wio **h dp ¥-^ 

a 2 V0 (T) ,oc J() ^ r) 



2tt J dp p i-2m(T)i ( 44 ) 

where Jo(?/) is a Bessel function of the first kind. As before, we use a ultraviolet cutoff 
A = ix j a to evaluate G(0): 

T^o (T) 

G(0) = W (45) 



We obtain finally 



gm cm) - 1 f( 2a Y mm r fa(r)l ^ |T », (lfi . 

G(r)-G(0)--n-j >, (46) 



which for small t]o(T) can be rewritten as 



G(r) - G(0) ' 



4vrr/ (T) 



r / 



(47) 



and in the ?7o(^) ~" > limit Eq. (j4T?j) is obviously reproduced. 

For arbitrary values of < f]o{T) < 1 it is not obvious to see how the KT theory is 
recovered when the anomalous phonon spectrum is taken into account. Recall that the 
logarithmic behavior in Eq. (J4U|) is crucial in the KT argument in the presence of vortices. 
Indeed, a simple scaling argument with free energy of the vortices combined with the results 
of spin- wave theory allows us to determine the value of t)o{T) at the critical temperature 
T c .— A more elaborate argumen1j^2i using the RG shows that the stiffeness is renormalized 
and i]o{T) becomes rj(T) = Tm 2 /2np s {T). However, the value of rj(T) at T c is the same 
as the value of 1] (T) at T c . Such an RG analysis led to the cellebrated prediction of a 
universal jump of p s (T) as T c is approached from below.— At first sight we may think that 
an anomalous phonon spectrum would disrupt the whole argumentation due to the form of 
the corresponding non-local hydrodynamics (|3ip. In the following we will show that this is 
not the case. To this end it is necessary to proceed in two steps. First, we make an analysis 
about the validity of the hydrodynamic description given by flHlj) . which will provide us with 
a lower bound for t)q(T). Second, we consider a careful analysis of the vortex field theory 
associated to this problem. We will see that the anomalous contribution to the vortex field 
theory becomes irrelevant at large distances and that, effectively, the same vortex field theory 
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as in the KT case holds. This result will in turn provide us with a physical interpretation 
of the lower bound for r] (T). 

Let us then start by discussing in more detail the effective Hamiltonian (J31j) . In order to 
define the mass density (132)) we have to be able to Fourier transform M.{p) = nm{pa)~ 2 ' m<KT \ 
Such a Fourier transformation is only possible if 1/4 < i] (T) < 1. To see this we consider 
the integral 

/"2-7T />oo gipr cos </) 

-^( r ) = /<-, \? / / dp 



2tt y„ ' p>™m-i 



2 n r 2[l- Vo (T)] J o y y 2 V0 (T)-l- 

The power of r in the second line of Eq. ()48|) must be positive in order to make the above 
result well defined at large distances. This gives us once more the upper bound for r] (T), 
i.e., T]o(T) < 1. The lower bound for f]o{T) follows from the asymptotic behavior of the 
Bessel function for y large, 

Mv) ~ \ — cos(3/-tt/4). (49) 
V 

Thus, the integrand in Eq. (|4l?j) behaves for large y like ~ i/y 2j M T )- 1 / 2 anc j ft follows 
that in order to avoid a power- like ultraviolet divergence we must have 2i]q(T) — 1/2 > 
0, which leads to f]o{T) > 1/4. Therefore, in order to have a well defined mass density 
A4(r — r') the inequality 1/4 < f]o{T) < 1 has to be fulfilled. Note that the upper bound is 
associated to a large distance (infrared) divergence while the lower bound is to short distance 
(ultraviolet) divergence. We have already seen that saturation of the upper bound leads to 
the determination of the temperature T*. On the other hand, we will give arguments below 
showing that saturation of the lower bound determines the actual critical temperature of 
the system. 

At higher temperatures, near the phase transition, we are allowed to keep only the zero 
Matsubara mode, such that the field theory becomes two-dimensional. The vortices are in- 
troduced in the standard way,^ and a duality transformation gives the following generalized 
sine-Gordon action for the dual field theory of vortices: 

" T 



Sduai = / d z r / d 2 r' 



In 2 



r x (r — r')d r (f(r) ■ d r icp(r') — z 5 2 (v') cos^(r) 
13 



(50) 



where T(r — r') = mn5 2 (r — r') + A4(r — r') and z is the fugacity of the gas of point vor- 
tices. A derivation of the action (|5Uj) is given in Appendix B. Nonlocal gradient terms in 
sine- Gordon theory have recently been discussed in a different context , 22 > 23 where a renor- 
malization group (RG) analysis shows that a local gradient term (Vy) 2 with a positive 
coefficient is generated.— The nonlocal contributions can be rewritten as a sum of higher- 
gradient terms which are all irrelevant at large distances. Let us state this in more simple 
terms. In momentum space the non-local gradient has the form 

a 1 I d P p 2[i +% (T)]^ (pM _ p) . (51) 



87T 2 nm J (2ir) 2 

Thus, as a local gradient term is generated, we have that for small p (large distances) the 
p 2 -term dominates over the p 2 ^ 1+r io( T )} one anc [ the anomalous contribution can be neglected. 
This argument shows that the dual theory (|5U|) actually has a KT transition, since the 
fluctuation generated p 2 -term is dominant in the infrared, leading effectively to a sine- 
Gordon theory of the usual type. Thus, our theory will ultimately be in agreement with 
the analysis of Ref. [l6| where no anomalous dimension arising from RPA corrections is 
considered. 

In the low-temperature phase, where < i]q(T) < 1/4, the hydrodynamic description via 
the effective Hamiltonian (jUlj) breaks down, since the inverse Fourier transform of A4(p) is 
no longer defined. This regime is equivalent to the one in the dual theory where the non-local 
gradient term becomes irrelevant. This means that the effective Hamiltonian of the classical 
theory in this range of temperatures is actually given by ([3D)). In other words, when the 
anomalous sine-Gordon theory becomes at large distances effectively the usual sine-Gordon 
theory (remember that a local gradient term is generated by fluctuations), we can dualize 
it back to obtain the effective original theory as given by Eq. (jHOJ)- 

In a KT transition the critical temperature T c is determined from the equation f](T c ) = 
1/4, or equivalently, T c = 7rp s (T c )/(2m 2 ).- Note that f]o(T) corresponds to a low-temperature 
approximation to i](T). Thus, we can interpret that saturating the lower bound in the 
spin- wave inequality 1/4 < r] (T) < 1 leads to an approximate value of the actual critical 
temperature for the KT transition. Therefore, we have that 

,52) 

We want to emphasize that the agreement of the value of i]o(T) at the lower bound of the 
inequality 1/4 < r] Q (T) < 1 with the value r](T c ) = 1/4 obtained from the KT theory is not 
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a simple coincidence. It follows from the fact that the temperature at which f]o(T) = 1/4 
corresponds to the onset of the regime where the non-local term of the generalized sine- 
Gordon theory becomes irrelevant. 

We have obtained tjq(T) as a correction to the power of the excitation spectrum at finite 
temperature by using a classical approximation to evaluate the density correlation function. 
Thus, we have simply accounted for a classical effect in a quantum calculation, i.e., our 
calculation of the spectrum at finite temperature is actually a semi-classical one. In the 
KT theory the value rj(T) = 1/4 is follows due to a conspiracy between spin-wave theory 
and the statistical mechanics of vortices.— In our case, a similar result obtains: our bound 
for f]o{T) is derived by analysing the non-local spin- wave theory, while the conclusion that 
r)o(T) = 1/4 determines the critical temperature needs additional analysis involving the 
vortex field theory (|50p. A more technical explanation follows by recalling that density and 
longitudinal phase fluctuations are related through the Ward identities.— More precisely, in 
Ref. |2J] it is shown that the density correlation function is related through an exact identity 
to the longitudinal component of the current correlation function, which corresponds to the 
response of the system to the longitudinal phase fluctuations or, in other words, spin-waves. 
Since 1] (T) is determined by the density correlation function, by taking the classical limit of 
this identity we obtain a relation between r]o(T) and the spin- wave response. The vortices, 
on the other hand, are related to the transverse phase fluctuations. It can be shown^ 
that the transverse component of the current correlation function decouples from all other 
correlation functions. This result holds also in the classical limit and that is the reason why 
the classical statistical mechanics of vortices alone can determine the phase transition. 

From the lower and upper bound to f]o{T) we obtain that the actual T c obtained from 
the KT theory is smaller than T*. Indeed, we have that 

rp = nPs(T c ) < 7rp B (0) < 2np s (0) = ^ 
2m 2 2m 2 m 2 

Since T c < T*, it seems to exist a region T c < T < T* where the system is not a superfluid, 
while not being a normal fluid either, due to the anomalous exponent. This region would be 
actually very small: from the inequality (}23|) we see that a more correct statement would be 
to claim the existence of a temperature region T c < T <C T* where an anomalous normal state 
occurs. This state can be thought as a quasi-condensate without phase coherence. In this 
respect, it bears some resemblance with the pseudogap state in high-T c superconductors,— in 
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which the phase fluctuations above the superconducting critical temperature play a similar 
role. 25,26 In the pseudogap phase the spectrum is not like the one of a normal metal, although 
the system is not in a superconducting phase. However, we should warn the reader that the 
actual physics of the pseudogap state is likely much more complex and that such an analogy 
must be considered with utmost caution. 

V. CONCLUSION 

In this paper we have shown that temperature effects induce roton-like excitations in a 
dilute two-dimensional Bose gas. From this result it followed that for nonzero low tempera- 
tures the spectrum is not of the phonon type and has an anomalous scaling with temperature 
dependent exponent. Thus, we have obtained in the quantum regime a situation which is 
reminiscent of the Kosterlitz-Thouless transition, namely, a continuously varying exponent. 

It would be interesting to extend this analysis to strongly-coupled two-dimensional Bose 
systems such as films of 4 He, where experimental data have so far been fitted by a superfluid 
density whose vortex-independent background contribution is calculated from the Bogoli- 
ubov spectrum.— It was found that the phonon contribution alone with its T 3 behavior 
(recall that for a phonon spectrum we find in general a contribution behaving like T d+1 in d 
dimensions) is not consistent with the temperature dependence of the data, thus calling for 
an improvement of the theory. In order to apply our approach to 4 He films we must derive 
the t- matrix for the actual atomic interaction in helium beyond the dilute limit. Prelimi- 
nary work has been done some time ago,— and phonon as well as roton excitations of the 
spectrum were obtained. However, these works concentrate only on the low-temperature 
properties. We believe that thermally induced roton-like excitations should also occur in 
this case, though in a more complicate manner. Certainly, such strong-coupling problems 
require more powerful calculation methods, for example, field-theoretic variational pertur- 
bation theory as developed in Ref. 12, which has led to the most accurate predictions of 
critical exponents so far.— 
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APPENDIX A: CLASSICAL LIMIT OF THE POLARIZATION BUBBLE 
Eq. (|11|) can be rewritten as 

U(iu, p) = 4m 2 / -^-n B (■£-)( \ \ 1 . (Al) 

v ^ ! J (2vr) d \2mJ \2miuo - p 2 - 2p • q 2miu + p 2 - 2p • q ) v ; 

In the classical approximation we write n^x) ~ T/x and the polarization bubble can be 
rewritten as 

U(icu, p) = 4m 2 T(/ + - (A2) 

where 

j = _ % f _A 1 1 (A3) 

The integrals I± can be evaluated using the Feynman parameters 

»oo poo M 
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l± = -I 



poo poo M n 

\ dXl / rf A 27 ^ 7 e- Al(2mw±ip2+2ip - q) e- A292 . (A4) 
Jo Jo ( 27r ) 



After evaluating the Gaussian integral over q we obtain 

/ ,/\, / ,/.V, I — 

(27T) d 



poo poo / \ ' 



T { ^T(d/2 - l)r(8 - <i) P « (l =F . (A5) 



2 d n d/2 V / /V V P 2 / 

Substituting the above expression back into (jA2|) we obtain Eq. (|T5j). 

APPENDIX B: DERIVATION OF THE ANOMALOUS SINE-GORDON AC- 
TION 

At sufficiently high temperatures the Matsubara time dependence of the phase 9(t, r) can 
be neglected and the effective action (|2*H|) can be written simply as 



Sett ~ ^ / rfV / rf 2 r'Ai (r - r') v s (r) ■ v s (r'). ( Bl ) 
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Although the duality transformation can also be performed in the continuum, a more tech- 
nically correct analysis is obtained in the lattice formalism. 4 The lattice version of the above 
action suitable for a duality transformation is given by the Villain form 



J)i 



(B2) 



where we have set the lattice spacing to unit and the components of are understood 
as lattice derivatives. The field rij is an integer field defined on the lattice. The partition 
function is then given by 

Z = T£Ut^MS L ). (B3) 

{n;} n i 

The first step in the duality transformation is the introduction of an auxiliary field through 
a Gaussian completion, i.e., 



exp 



2m 2 T 



oo 



-v / Y\ dbkp, exp 

- 50 fe lAt 



m 2 T 



E HM-%b 3 - i E V (VOj ~ 2™,- 



(B4) 



Next we apply the Poisson formula 

oo /»oo 

E / dxf(x)e i2nnx 



E /h> 



(B5) 



m=— oo 



in order to convert the integral over the real auxiliary b; field into a sum over an integer 
field Nj. After this manipulation the periodic field 6i can be easily integrated out after a 
summation by parts. This leads to the constraint V • Nj = 0. Up to unimportant overall 
factors the partition function becomes 

m 2 T 



Z = E <W 4 ,oexp 



(B6) 



In two dimensions the constraint is solved through 



Ni 



ifj, 



(B7) 



The resulting partition function will be the one of a neutral Coulomb gas in the lattice. By 
applying the Poisson formula once more, we obtain 

m 2 T 



E / n rfv9feexp 

{ Sl } J -°° k 



E v»<Pi(M ^ijVpipj - i E 27rs m 



i-.i 



(B8) 
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The integer fields Sj represent the point vortices of the theory. We can now introduce 
the vortex fugacity z to build a grand-canonical ensemble of vortices.— The most relevant 
vortex configurations correspond to Sj = ±1. After the rescaling ipj — > ifj/2ir and taking 
the continuum limit of the grand-canonical ensemble theory, we obtain Eq. (J5(Jj) . 
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